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Abstract: We reformulate and extend our recently introduced quantum kinetic theory for 
interacting fermion and scalar fields. Our formalism is based on the coherent quasiparticle 
approximation (cQPA) where nonlocal coherence information is encoded in new spectral 
solutions at off-shell momenta. We derive explicit forms for the cQPA propagators in the 
homogeneous background and show that the collision integrals involving the new coherence 
propagators need to be resummed to all orders in gradient expansion. We perform this 
resummation and derive generalized momentum space Feynman rules including coherent 
propagators and modified vertex rules for a Yukawa interaction. As a result we are able to 
set up self- consistent quantum Boltzmann equations for both fermion and scalar fields. We 
present several examples of diagrammatic calculations and numerical applications including 
a simple toy model for coherent baryogenesis. 
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1. Introduction 

Many interesting problems in particle physics and cosmology involve coherent, relativis- 
tic quantum systems in out-of-equilibrium conditions, examples including the electroweak 
baryogenesis [1-5], leptogenesis [6-9] and the out-of equilibrium particle production [10]. 
While it may be straightforward to give a mathematical formulation for a given problem, 
for example using the Schwinger-Keldysh method, it is often hard to find an approxima- 
tion scheme that is accurate, but simple enough to be used in practical applications. We 
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have recently developed a new transport theory for relativistic quantum systems including 
nonlocal coherence [11—13]. Our formalism is based on the observation [11] that in systems 
with certain space-time symmetries the nonlocal coherence information is encoded in new 
spectral shell solutions for the dynamical 2-point functions at particular off-shell momenta. 
The extended spectral structure of the two-point functions, the coherent quasiparticle ap- 
proximation (cQPA) [14,15], was first discovered in [11] and a complete transport theory 
with collision terms for fermions was presented in [12] and for scalar fields in [13]. 

In refs. [11, 12] the fermionic transport theory was formulated using a convenient, but 
somewhat less familiar spin-projected 2-component notation. Later, in refs. [14, 15] it was 
noted that the collision integrals involving the coherence functions need to be resummed to 
all orders of the gradient expansion. In this paper we reformulate our formalism for fermions 
using a more familiar 4-component Dirac notation, accompanied by a much more convenient 
parametrization of the coherence-shell solutions. We will also show that the resummation 
issue is much more generic than what was considered in [14,15]: a resummation is necessary 
for all external and internal lines involving coherence propagators in the diagrammatic 
expansion of all collision integrals appearing in the transport equations. Here we will show 
how this resummation can be done such that all effects are absorbed into a redefinition 
of the effective coherence propagators. Using these results we develop an extension of the 
momentum space Feynman rules for diagrammatic evaluation of self-energies and eventually 
collision terms including quantum coherence effects. We will give a full discussion of single 
flavour fermion and real scalar fields, interacting via a Yukawa type of interaction. In 
addition to our extension of the usual Feynman rules to the coherent case, the main results 
of this paper include a derivation of simple, generic quantum Boltzmann equations (qBE) 
for fermions and scalars. Fermionic equations couple the usual mass-shell distribution 
functions to new coherence functions with compact expressions for the collision integrals 
as traces over Dirac- matrix projections of generic self-energy functions. The scalar field 
qBE's on the other hand are cast into a form of a finite set of moment equations, again 
with a well defined collision integral that can be evaluated diagrammatically by use of our 
cQPA Feynman rules. 

We use our generalized Feynman rules to compute several perturbative self-energy dia- 
grams in the Yukawa theory, including the full coherent fermion and scalar field self-energies 
at the one loop level. We also study a 2-loop example which bears resemblance to a more 
involved calculation with flavour mixing fields [16], necessary for a resonant leptogenesis 
problem. These examples illustrate the ease of use of our formalism for diagrammatic eval- 
uation of collision integrals involving the nonlocal quantum coherence effects. We also com- 
pute explicit expressions for the collision integrals following from the Yukawa interaction at 
the one- loop level and interpret these results in terms of the mass-shell and coherence-shell 
contributions in scattering (decay and inverse decay) processes. As a numerical applica- 
tion of our formalism, we solve the cQPA phase space distribution functions for a fermion 
with a complex, homogeneous time varying mass parameter evolving smoothly from one 
phase to another, in presence of decohering interactions. These results illustrate the ef- 
fect of collisions on the production and evolution of the coherence and particle numbers 
in the presence of interactions. We also build a toy model for coherent baryogenesis (see 
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eg. ref. [17]) by assuming that the coherently evolving fermion ip decays into the ordinary 
standard model quarks. We compute the temporal evolution of the seed quark asymmetry 
and the ensuing baryon asymmetry and show that such a scenario, when embedded to a 
more complete particle physics model, could be an attractive mechanism for baryogenesis. 
In this paper we only consider the spatially homogeneous and isotropic problems. 

This paper is organized as follows. In section 2 we briefly review the derivation of 
the fermionic cQPA shell structure, and the equations of motion that follow from using 
the cQPA structure as an ansatz in the full dynamical Kadanoff-Baym equations. We 
also introduce a more convenient parametrization of the coherence propagators using the 
familiar 4-dimensional Dirac notation with projection operators. In section 2.4 we present 
the resummation method giving rise to a consistent expansion of the fermion collision 
integrals. In section 3 we repeat the analysis for scalar fields, and in section 4 we derive 
the generalized Feynman rules for the computation of the resummed self-energies for a 
Yukawa-theory. In section 5 we compute examples of the self-energy diagrams and collision 
integrals and in section 6 we apply the formalism to a toy-model of a coherent baryogenesis. 
Finally, section 7 contains our conclusions and outlook. 

2. cQPA for Fermions 

We will be using the Schwinger-Keldysh approach [18] to non-equilibrium quantum field 
theory (for a review see ref. [19]). The basic objects of interest for us are the dynamical 
Wightman functions 1 S < ' > : 

iS < {u,v) = {i)(v)^(u)) and iS>(u,v) = (ip(u)ij>(v)) , (2.1) 

where (...) = Tr{/5...} denote the expectation values w.r.t. some unknown quantum 
density operator p. The time-ordered (Feynman) and anti-time-ordered (anti-Feynman) 
2-point functions are then related to the Wightman functions by 

S F (u, v) = 9(u - v )S > (u, v) - 9(v - u )S < (u, v) 

S p {u, v) = 9{v - u )S>(u, v) - 9(u - v )S < (u, v) . (2.2) 

The basic idea in the SK-approach is to write down a (Schwinger-Dyson) hierarchy of 
equations for all Greens functions in the theory, thus bypassing the need to define p at all. 
The main approximation is a truncation of this hierarchy to a closed set of equations for a 
restricted set of correlators. In practice the truncation is often done already at the level of 
the 2-point functions, assuming that the self-energy functions in the SK-equations can be 
evaluated perturbatively. In this case only the two Wightman functions (2.1) are needed 
get the closure 2 . It is convenient to separate the internal and external scales by performing 
Wigner transformations to the 2-point functions: 

S(k,x) = J d 4 re ik - r S(x + §,x- §) , (2.3) 

x Note that our sign-convention for S < is opposite to the usual one, see eg. ref. [5]. 

2 Given the Wightman functions, the spectral function is found from A = |(5' > + S < ). After this the 
spectral relation Sh(u,v) = — isgn(u° — v°)A(u,v) fixes the retarded and advanced correlators S r , a = 
Sh ± iA, which completes the maximal set of independent 2-point functions. 
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where x = (u + v)/2 is the average coordinate, and k is the internal momentum conjugate 
to the relative coordinate r = u — v. It is straightforward to show [5, 12] that in the mixed 
representation the correlators S < ' > satisfy the Kadanoff-Baym equations [20]: 

(ft+ l -p x -mo - im^SO - e-^{^ H }{S<'>} - e - l< >{Z<'>}{S H } = ±C colh (2.4) 
where "+" refers to "<" and "— " to ">" . The collision term is given by 

Cecil = ({Z>}{S<} - {£<}{£>}) , (2.5) 

where the 0-operator is the following generalization of the Poisson brackets: 

0{f}{ 9 } = \ [d x f • d k g - d k f • d x g] . (2.6) 

We are interested in problems where a complex mass function can be spatially or temporally 
varying m = mn(x) + iraj{x), whereby the nontrivial mass operators mo and appear 
in Eq. (2.4). Their action on ^-functions is defined as: 

m 0i5 S(k, x) = m RJ (x)e-^- 9 is(k, x) . (2.7) 

The self-energies £ are some complicated functional of the correlators S < ' > , and they 
can be computed using some explicit truncation scheme, such as the 2PI effective action 
in loop expansion [19,21]. From Eq. (2.4) we immediately see that the spectral function 
A = ^{S > + S' < ) obeys an equation identical to Eq. (2.4), but with a vanishing collision 
term. In addition the spectral function needs to satisfy the sum-rule: 

J ^U(M)7° = 1, (2-8) 

which follows for example from the equal time anti-commutation relations of the fermionic 
fields. We will later see that this constraint is strong enough to completely fix the spectral 
function in our approximation scheme. Equations (2.4) for S < and A, together with the 
sum-rule (2.8) form a complete set of equations when the interactions (a scheme to compute 
S) and the mass profiles are specified. In practice, these equations are too hard to be 
solved in their full generality, and several approximations are needed to find a solvable set 
of equations. 

2.1 Approximations 

The usual kinetic approach to non-equilibrium quantum field theory [5, 19, 22] considers 
weak interactions and slowly varying classical backgrounds. In this case the excitation 
spectrum of the system can be reasonably approximated by a singular shell structure and 
the KB-cquations (2.4) reduce to Boltzmann equations for the on-shell excitations called 
quasiparticles. Moreover, the standard approach typically makes the assumption that the 
system is, as a consequence of being close to a thermal equilibrium, nearly translation 
invariant. This implies that the phase space singularity structure is also essentially the 
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same as in thermal equilibrium. The key idea leading to the coherent quasiparticle approx- 
imation [11-15] scheme is to relinquish the assumption of nearly translation invariance. 
The ensuing new solutions which describe the nonlocal quantum coherence are typically 
oscillatory in the space and time coordinates in quantum scales x ~ fe -1 , which explicitly 
violates the translation invariance. To find the cQPA we make the following approximations 
to Eqs. (2.4): 

1) We neglect the terms oc Sh on the LHS of the KB-equations (2.4). 

2) We neglect the terms oc T,h on the LHS of the KB-equations (2.4). 

3) We consider the spatially homogeneous and isotropic case with V5' < = 0. 

The first approximation above is made also in the standard kinetic theory and it should 
be reasonable at least in the limit of weak interactions. It has been shown that for a 
scalar field close to thermal equilibrium the neglected terms are of higher order in the 
interaction rate T than the dominant contribution from the collision term [5]. From the 
practical point of view, these terms must be neglected in order to find a singular phase space 
structure for the spectral function, even in the usual Boltzmann approach. The second 
approximation is made only for simplicity; the E# corrections would merely give rise to 
modified dispersion relations for the quasiparticles, but not change the generic structure 
of the theory. These terms could be incorporated later on. The third assumption on the 
other hand is crucial. The simple spectral structures for the coherence solutions, which 
is the heart of the cQPA, arise only in systems with particular space-time symmetries, 
including spatially homogeneous and isotropic, and stationary planar symmetric cases [11]. 
This paper is devoted to the former of these symmetries. 

With the approximations 1-3 the KB-equations (2.4) can be reduced and decomposed 
into the following Hermitian (H) and anti-Hermitian (AH) parts: 

(H) : 2k S< = HS< + S<& + ij°(C co n - C^j (2.9) 
(AH) : id t S< = HS< - S<H* + ij°(C coll + cl n )j° , (2.10) 

where we defined a hermitian Wightman function S < = iS < "j°, and the operator 

H = k • a + 7°m + ij°j 5 m 5 (2.11) 

can be interpreted as a local free field Hamiltonian. Because of the spatial homogeneity 
the mass functions rriRj{t) depend now only on time. Note that in the collisionless limit 
(C C oll ~~ * 0) and in the zeroth order in gradients (rriQ^ — > mm) only the AH-equation 
contains an explicit time derivative of S < , while the hermitian equation becomes a purely 
algebraic matrix equation. For this reason Eq. (2.10) is called a "kinetic equation", and 
it describes the time evolution of the Wightman function S < . The role of the hermitian 
"constraint equation" (2.9) on the other hand is to restrict the phase space structure of 
S<. 

Even after approximations 1-3 the resulting equations (2.9-2.10) remain formidable to 
solve. Note in particular that they still involve the mass operators containing terms which 
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are of arbitrarily high order in gradients. In order to find a tractable set of equations we 
will make the following two additional approximations: 

4) We will expand the non-dynamical KB-equation (2.9) to the zeroth order in the 
scattering width V = |(S > + S < ) and the mass gradients dtm. This gives rise to the 
singular cQPA phase space structure. 

5) We insert the singular cQPA structure as an ansatzto the dynamical KB-equation (2.10), 
and expand also this equation to the lowest nontrivial order in T and dtm. 

As we shall find out in the subsection 2.2 below, the cQPA phase space structure contains 
two independent mass-shell distribution functions f m h± and two coherence distribution 
functions f c h± for each helicity h. When this shell structure is fed into the dynamical 
equation (2.10) and the equation is integrated over ko, it reduces to the sought after 
quantum Boltzmann equations for the on-shell functions f a ■ We shall derive these equations 
in the subsection 2.3. However, because the coherence-shell functions are rapidly oscillating 
in time, a special care is needed when expanding the collision terms in these equations. A 
novel resummation scheme that gives a consistent expansion of the collision integrals to 
leading order in T and dtm will be introduced in section 2.4. 

2.2 cQPA phase space structure 

The most general fermionic 2-point function S(k,t) consistent with spatial homogeneity 
and isotropy consists of the terms with products of 7 , k • 7 and 7 s only, and thus spans 8 
of 16 independent components of the full Dirac algebra. Furthermore, we observe that the 
helicity operator: h = (k ■ 7°77 5 )/|k|, commutes with the Hamiltonian (2.11), implying 
that different helicity projections of S < do not mix in a collisionless theory, i.e. helicity is 
a good quantum number. For this reason, it is convenient to decompose the 8 independent 
components of the 2-point function S < (k, t) in terms of helicity projectors: Ph = \{\ + hh), 
h = ±1 to get: 

iS < (k,t) = \Y1 p hh°9h0 + hk- 19h3 + g hl + i^g h2 ] , (2.12) 
h=±i 

where gh a (k, t) are the 8 real scalar components. Using the the chiral (Weyl) representation 
of the Dirac algebra: 7 = a 1 (g> 1, 7 = ia 2 0a and 7 5 = — a 3 ®l, the helicity block-diagonal 
Wightman function in Eq. (2.12) can be written as 

S<(k,t) = £ gZ(k,t)®±(l + hl-v), (2.13) 
h=±l 

where the Bloch representation of chiral part is just g"^ = \ (g^ Q + g^cr 1 ) ■ 

We now proceed to step 4 in our approximation scheme and analyze the constraint 
equation (2.9) in the zeroth order in T ~ T, <,:> (neglecting the collision term) and dtm. In- 
serting this decomposition into the constraint equation (2.9) reduces it to two independent 
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homogeneous matrix equations for the helicities h = ±1 (note that the index ordering is 
here defined as a = 0, 3, 1, 2): 



Equation (2.14) may have nontrivial solutions only if 
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(2.14) 



det(B h ) = (k 2 - \m\ 2 ) k 2 = 0. 



(2.15) 



This constraint gives rise to a singular shell structure, because the solutions need to be 
proportional to either S(k 2 — k 2 — \m\ 2 ) or 5(ko). The full matrix structure of the cor- 
responding singular solutions can then be worked out by setting ko 7^ and ko = in 
Eq. (2.14), respectively, so that the combined solution is eventually found to be [11,12]: 



is<(k,t) = j2 p h 



■t 



27rsgn(/e )(^+ m R - il 5 mi)f< hsgn{ko) 5(k 2 - \m\ 2 ) 



+7T k ■ 7 



J hi Jh2 



V Ikl 



fhl ~ ^fh2 W 



(2.16) 



Here the first line is the standard one-particle mass-shell solution with dispersion relations 
ko = it^k = ±(k 2 + |m| 2 ) 1 / 2 , while the second line is a new coherence solution living at shell 
ko = 0. The real mass-shell distribution functions /^^(k, t) can be related to the phase 
space densities for particles and antiparticles via the Feynman-Stiickelberg interpretation 3 : 



ikfc(t) = fmh+(k,t) n^ h {t) = 1 - /< h _(-k,t) 



(2.17) 



and the new (real) ko = 0-shell distribution functions f?i 2 (k, t) describe nonlocal quantum 
coherence between these modes. The representation (2.16) for S K was derived in rcfs. [11, 
12]. For our purposes here it is actually more convenient to reparametrize the ko = 
coherence-shell functions as follows: 



f * h± ~ 2(w 2 



m 2 R ) 



(2.18) 



Using these new functions the Wightman function Eq. (2.16) becomes 
iS<(k,t) = 2irY^±-}-P h (k' ± + m R -i 7 5 mi) 



h± 



/^%T^)+(7°T^)4 



S(k ) 



(2.19) 



where we denote k± = (±Wk,k). The complete singular phase space structure (2.19) is the 
coherent quasiparticle approximation for S < . 



3 We associate the negative energy states with the CP- rather than C conjugates to the positive energy 

— / \ op 

states: n kh = (n kh ) . 
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A note on conventions. The shell functions appearing in Eq. (2.19) are labelled by three 
indices: j = m,c referring either to mass- or coherence-shell, h = ±1 referring to helicity, 
and ± referring to the positive and negative energy. These indices will be appearing, always 
in this order, in many other constructs throughout the paper. Sometimes one or the other 
of the indices is to be summed over, and in order to suppress the notational clutter we will 
always use an implied summing convention, where a non-displayed index is summed over. 
That is, we use for example 

E Qj*± = E = E Qi = Q > ( 2 - 2 °) 

jh± jh j 

where Q can be for example S <:> , f K,> or any other construct defined using them. 

Note that the time dependence of the on-shell functions f a is not defined yet; in the 
next section we will use the cQPA structure as an ansatz in the dynamical equation (2.10) 
to derive the generalized quantum Boltzmann equations for f a . First however, we have to 
determine the cQPA propagators also for the function S > , the spectral function and the 
pole propagators. 

2.2.1 Pole functions and thermal limit 

In the collisionless limit S > and the spectral function A = ^(S > + S < ) satisfy the same 
constraint and kinetic equations as S < . The spectral solution must then be identical to 
Eq. (2.19) with four yet undefined on-shell functions f& for both helicities. However, the 
spectral function must in addition satisfy the sum rule (2.8), which completely fixes these 
functions: f^ h± = \ and f^ ± = [11]. The spectral function then reduces to the familiar 
form: 

A(k, t) = 7rsgn(/c )(^ + m R - ij 5 rrn)S(k 2 - \m\ 2 ) . (2.21) 

The retarded and advanced propagators are directly related to the spectral function: 
S r ,a(u,v) = z f2id(±(uo — vo))A(u, v), so that they are also given by their standard ex- 
pressions: 

W ) / + " 5 7' . ( 2 . 22 ) 
k z — \m\ A ± 2sgn(/eo)e 

Next, the relation A = |(5* > + 6' < ) together with Eqs. (2.19) and (2.21) completely specify 
S* > ; it has an expression identical to Eq. (2.21), where the on shell functions are given by 

/^± = l-/^± ^d f> ± = -f< ± . (2.23) 

Thus, the four on-shell functions f£ are the only dynamical variables in the cQPA; all 
other 2-point functions can be expressed in terms of these. For this reason we only need 
to consider the dynamical equations for S < and the on-shell functions /<. 

Let us stress that the new coherence solutions appear only in the dynamical functions 
5 <5> in the cQPA, but not in the spectral function A or in the pole propagators SV,a- This 
is quite natural, as it merely expresses the fact that there cannot be coherence without the 
interfering mass-shell excitations; coherence is a purely dynamical phenomenon. Of course 
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the Feynman and anti-Feynman propagators do have well defined coherence contributions 
that can be computed from their respective definitions: 

Sp = S r — 

Sp = -S a - S< . (2.24) 

To conclude this section we note that in thermal equilibrium limit the full translation invari- 
ance kills the coherence-shell contributions also from the Wightman functions S ,<,> [11]. 
Moreover, the thermal density operator implies Kubo-Martin-Schwinger (KMS) relation 
[23]: 5>(fco) = e^ fc °5<(jfco), which fixes the remaining mass-shell functions to Fermi-Dirac 
distribution: f^ h± — > / eq (±Wk) = l/(e ±/3tUk + 1). This reduces S K to the standard thermal 
form [24]: 

iS< = 2vr sgn(A: )(^ + m R - i 7 5 mj) f cq (k ) 5{k 2 - \m\ 2 ) . (2.25) 

The thermal limits of the other standard propagators can then be found by using the 
general relations given above. 

2.3 Quantum Boltzmann equations 

As the final step in our approximation scheme, we insert the spectral cQPA-solution (2.19) 
as an ansatz into the kinetic equation (2.10). This will lead to a closed set of equations of 
motion for the dynamical variables f^ a . Let us first define a useful shorthand notation: 

iS< ^ 2vr ]T (S< h ±6(k T u*) + S< ± S(k )) , (2.26) 

h± 

where the spin-matrix functions Sj} l ± are easily read off from Eq. (2.19). Inserting this 
form to the kinetic equation (2.10) and integrating over ko we find the master equation: 

d t S< = -i[H, S<] + 7°(C coll + Cj oll )7 , (2.27) 

where we denoted (...) = J ^ (. . .) and, in accordance with the notation (2.26): 

(S<)=S<=J2(S< h± +S< ± ), (2.28) 

h± 

where = S^~ h ^°. Note that all fco-derivatives in the operator H have vanished in the 
integration process as total derivatives, reducing H to the usual local Hamiltonian: 

H = k • a + 7°m R + i^^rm . (2.29) 

The master equation (2.27) is in fact the most compact form of the quantum transport 
equations in the cQPA scheme. One can see that the lowest moment of the correlation 
function, (S < ) = 5 < , corresponds to a density matrix in the Dirac indices, and what 
remains to be done is to express the collision term in Eq. (2.27) in terms of the degrees 
of freedom included in S < . This is of course not possible in general, but it can be done 
in the context of the spectral cQPA ansatz. In the case of fermions it is convenient to 
go further than just computing the collision matrices however, and rewrite Eq. (2.27) as 
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a set of scalar Boltzmann equations for the on-shell functions f a , in exact analogy to the 
ordinary Boltzmann equation. 

We begin by writing the integrated correlation function S < as (we will drop the su- 
perscript < on f a when there is no risk of confusion) : 



S < — PhPk± fmh± + 7tk±fch± 



h± 



where 



1 



1 ± — and 7£ k± = 7 U =F 



_ 0^ m R 



(2.30) 



(2.31) 



Here P k -t are positive and negative energy projection operators. Indeed, using the the 
identity H 2 = lo^ one can easily show that HP\r± = ±u; k P k -|- and then 

Pk+Pk- = Pk-Pk+ = and Pl ± = P k± . (2.32) 

Moreover, the coherence matrix TZ^± and energy projectors P k -t obey the relations: 

P k± ft k± = ft kT P kT = P k±7 °^k T , (2.33) 

where the last equality follows from the anticommutator {H, 7 } = 2m#. The matrix lZk± 
is not a projector, but combined with the energy projector it acts as an effective projector 
onto the coherence solution under the trace operation, such that the on-shell functions f a 
can be solved from Eq. (2.30) tracing over Dirac algebra: 



mh± 



Tr 



Pk±Ph<S < 



fch± = CkTr 



K k ±P k ±P h S < 



(2.34) 



where £ k = ^ k /(w k — rri^). The equations of motion for the on-shell functions fh a can 
now be obtained from the master equation (2.27) by forming the appropriate traces, or 
alternatively by taking time-derivatives of Eqs. (2.34) and using the kinetic equation (2.27) 
for dt<S < , the identities (2.32-2.33) and finally computing the traces over Dirac algebra. 
After some algebra one finds the sought after quantum Boltzmann equations: 



9tfmh± = ±2^W/ C ' l + + ^kh-fch-) + C m h±[f a ] 



9tfch± = T^kfch± + £k$kta — ~fch± - 7>Afmh + C c h±[f a ] , 

L ^k 2 J 



1 



(2.35) 
(2.36) 



where Af m h = f m h+ — fmh- and all dependence on the mass gradients is encoded in 



w (m R \i \\t\m\ 
\ w k / U) v 



After using cyclicity of trace the collision integrals Cjh± can be written as: 



Cmh±[fa] = Tr 7°(C coll + C c t oll )7° J Pk±P/ l 



Cch±[fa] = CkTr 7 (C coll + C c t oll ) 7 °^ k± P k± P /l 



(2.37) 

(2.38) 
(2.39) 
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Eqs. (2.35-2.36) together with the collision integrals (2.38-2.39) form a closed set of equa- 
tions of motion for the on-shell functions f a , once the self-energy functionals T, <:> ap- 
pearing in the collision integrals are specified. Beyond the introduction of the convenient 
4-dimensional notation, the main contribution of this paper is to show how the collision 
terms (2.38-2.39) can be consistently computed to a given (lowest) order in gradients and 
coupling constants. As we shall see, this is a nontrivial task due to the need for resumma- 
tion of the coherence propagators, which eventually leads to a new set of momentum space 
Feynman rules. 

2.4 Resummed fermion collision term 

The basic quantity appearing in the projected collision terms (2.38-2.39) is the zeroth 
moment of the collision term (C co ii). We shall now carefully explain how this quantity can 
be expanded consistently to the leading order in the interaction width and the gradient of 
the mass function, denoted by O 1 = 0(T,dtin). Let us begin with a straightforward, but 
limited in validity, derivation using the mixed space representation: 

(Ccoii) = ^ J ^e^({^>(k,t)}{S<(k,t)}-{^<(k,t)}{S>(k,t)}). (2.40) 

Naively, one would simply truncate the 0-expansion to the lowest order, as was done in 
refs. [11-13]. However, it was later realized that because of their oscillatory behaviour 
the coherence-shell functions contribute to the leading order term at all orders of the 0- 
expansion [14, 15]. We can see this quite easily when we observe that to order O 1 the 
equations of motion (2.35-2.36) imply that dtf m h± = and dtf c h± = T 2iwk/c/i±- This 
holds true also for the full mass-shell and coherence-shell Wightman functions (in the 
distribution sense): 

d t S< ± = O 1 , d t S<± = t 2iu; k S< ± + O 1 , (2.41) 

where S < = Yl±(^m± + ^c±) m accordance with the notation (2.20). If we make the 
additional assumption that c^E ~ O 2 , it is easy to show that (— iO) p {S > (/c, t)}{S < (k, t)} = 
((±uj] i dk ) p T l > )S^± + O 2 when p > 1. That is, we can use the naive expansion for the mass- 
shell part, but for the coherence-shell part we have to replace the — iO-operator by ±Wk<9fc 
acting only on the self-energy function S > : 

(Cecil) « \ E/ ^ (S>5< ± + (e^ d ^>)S< ± -(>♦*<)), (2.42) 

where we have dropped all terms of order O 2 or higher. We can now see that the exponential 
gradient operator in Eq. (2.42) precisely translates the self-energy functions associated with 
the coherence solutions S c ± from the actual ko = shell to the corresponding mass shells. 
This feature will persist also in the more accurate treatment where the smallness of dt'E is 
not assumed. 
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2.4.1 Consistent expansion and E e g 

Equation (2.42) consistently accounts for the resummation of all zeroth order gradient 
corrections from an external coherence-shell propagator in the collision integral. This 
introduces a correction to the naive gradient expansion of Eq. (2.40). It is easy to see 
that the same issue is encountered with the coherence propagators that appear in the 
perturbative expansion of the self-energy X. Indeed, whenever the time derivative acts 
on an internal coherence propagator in S, it returns as a leading term that propagator 
multiplied by ±2iu q , where oj q is the energy associated with the propagator. Such terms do 
not contain any gradient suppression and the assumption dfT, ~ O 2 made below Eq. (2.41) 
fails. Thus, we need a resummation scheme also for the coherence contributions inside S, 
before a consistent gradient expansion for (C co ii) can be established. We shall now formulate 
a general method of expanding (C co ii) consistently to the leading order. To this end, it is 
convenient to rewrite Eq. (2.40) in terms of the (partially transformed) coordinate space 
functions: 

(Ccon) = \Jdw [Z > (t,w ,k)S < (w ,t,k)-Z < (t,w ,k)S > (wo,t,k)] . (2.43) 

Note that this expression has only one external time t; forcing equal times exactly corre- 
sponds to integrating over the energy in the mixed representation. Because we are only 
interested in the spatially homogeneous case, all gradients vanish in the expression (2.43). 
Clearly, for the equations of motion (2.27) to close, we need to be able to express (C co n) en- 
tirely in terms of functions S^(k, t). We begin by writing the propagators in the two-time 
representation: 

S(w ,w' ,k) = J d 3 (w- w')e- ik < w -™">S(w,w') 

= f p. e -^ -<) s{koj k) !^±^o ) . (2 . 44) 

J 2-7T 2 

Inserting the spectral solution (2.26) to the definition (2.44) we get the two-time represen- 
tation of our cQPA-Wightman functions (t = (uq + vq)/2): 

iS<>>(uo,vo,-k) = £ [e^^-^S^M) +S<? (M) 

Note that iS < (t, t,k) = 5 < (k, t). We clearly need to expand the functions S < ' > (wq, t, k) 
appearing explicitly in Eq. (2.43). However, the perturbative expressions for the self-energy 
functions S <,> (t, wo, k) in general involve further integrations over internal vertices w' Q , w'q, 
and so we will need to expand correlators S(wq, w' , k) with arbitrary time-coordinates wq 
and w ; with respect to the "external" time t = (uq + vq)/2. This is easily done by Taylor 
expanding correlators around t: 

S<>>(k , k, ^) = £ 1 _ t y d?S <>>(k , k, t) . (2.46) 

n=0 



(2.45) 
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We can determine the derivative terms in the above equation by recursively using the zeroth 
order equations (2.41) for n > 1: 

d?S< ± = O 1 , d?S< ± = (Ti2^) n S< ± + O 1 . (2.47) 

This implies that the Taylor expansion of the mass-shell functions is trivial to the lowest 
order in mass gradients, and that the coherence solutions get multiplied by a simple phase 
factor. Inserting these results to Eq. (2.44) one finds: 

iS*'* (w ,v/ ,k) « ^2 [e^^-Os^g (k,t) + e Ti ^ W0+ <^ (k, i)] , (2.48) 

where we dropped terms of order O l . As a trivial consistency check we see that this 
expansion reduces to Eq. (2.45) as wq — > uq and w' — > vq. Physically the propaga- 
tor (2.48) is an approximation which takes into account the rapid temporal variations 
due to the oscillations of the coherence functions, but neglects the corrections due to 
temporal variations in the background fields. In particular equation (2.48) implies that 
iS < ' > (wo, t, k) = exp(=FKJk(^o — t))S± > , so that applying it to Eq. (2.43), we find 

(Cecil) = - \ £ (E> (io*)^ - £< (iu*)^) , (2.49) 

where 

?S>(fc.M = /d^C->^(^k). (2.50, 

The result (2.49) is similar to the equation (2.42), although in the new derivation we 
did not assume the smallness of <9t£. This is as expected, as both equations account 
for the resummation of the external coherence propagator. However, the correct self- 
energy function S e g in Eq. (2.50) is not in general equivalent with the projected mixed 
representation function S appearing in Eq. (2.42). Yet it is interesting to see that in the 
end, after the resummations have been carried out, Eq. (2.49) can be written in an integral 
form formally equivalent to the naive lowest order gradient expansion of the equation (2.40): 

(Ccoii)^/ ^(S> ff (fc,t)5<(fc,t)-S< ff (fc,t)5>(fc,t)). (2.51) 

The crucial point is that the original cQPA S <,> -functions have been replaced by the 
effective ones: 

iS<f (k,t)^2nY, (k, t) + S<±> (k, t) ) 5(fc T «*) = iS<>> + iS<£ . (2.52) 
± 

The resummed propagators iS^f* will also appear in perturbative expansions of the effec- 
tive self-energies . In section 4 we will return to this issue and derive a set of extended 
momentum space Feynman rules that can be used to formulate perturbative expansions 
and compute arbitrary self-energy diagrams in the cQPA framework. 

At the first sight it might look surprising that the resummations effectively pushed 
the coherence shells back to the mass shells in the effective propagators appearing in the 
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collision term. On the second thought this is of course perfectly natural: despite the fact 
that the ko = solutions are real (they have been constructed explicitly in exactly solvable 
free theories elsewhere [13,25]), they still physically correspond to coherence information 
between the on-shell states. What we are observing is our theory being consistent with 
the fact that it would make no sense to have a mere coherence information colliding with 
something; to have an interaction there always need to be a colliding particle present. When 
the coherence contributions are pushed to mass-shells their role is sensible: to modulate 
the collision frequencies of the mass-shell states depending on their state of coherence. 

2.4.2 Collision integrals for the on-shell functions 

Given the generic result Eq. (2.49) we can write the collision terms in Eqs. (2.38-2.39) 
explicitly as traces over the effective self-energy functions (2.50) multiplied by simple prop- 
agator matrices. To this end it is convenient to introduce some further notation, rewriting 
the coefficient matrices S m h± and S c h± as: 

S< h± (k,t) = ^-/C m ^±(k)/ mft± (k), 

5 ^± (M) = ^cfciW^k) , (2.53) 

where: 

KmhiQi) = (2cu k )P k ±P h (k) 7 ° = ±(j& ± + m R - z 7 5 m J )P /l (k) (2.54) 

and 

IC ch± (k) = (2u, k )P k± ft k± P ft (k) 7 ° = -*-IC mh± (k)IC mhT (k) . (2.55) 

2w k 

To derive the last form for K c h± we used the identity (2.33). Note that the operators IC m h± 
and )C c h± have the natural interpretations in terms of the spinor products. First, it is easy 
to show that 4 

IC mh ± = u±(k,h)u±(k,h) , (2.56) 

where u±(k,h) are the free positive and negative energy spinors normalized such that 
v}±u± = 2w k . From Eq. (2.55) one then immediately finds that the coherence propagators 
are proportional to off-diagonal spinor products: 

K ch ± = A h ±u± (k, h)u T (k, h) , (2.57) 

where the normalization factor 



A h± = -L M *, h)u±(ki h) = h M ± ^ _ m\ ™ . ( 2 .58) 
2u; k w k V w k / mr 



This observation gives rise to a natural normalization for the coherence-shell functions 
fch±'i t ne canonically normalized coherence functions clearly are: 



4 If one interprets our negative energy states as positive energy antiparticles by replacing k — > — k for 
the negative energy solutions one finds the familiar result: J2h K-mh± = fc±m, where m = rriR — i^°mi. 
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Figure 1: Graphical representation of Dirac structures of different on-shell propagators. Three 
momentum and fermion number flow from left to right. Arrows in fermion lines represent the 
fermion number flow and double arrows the direction of the energy flow. 



f, 



ch± 



Ah±fch± 



(2.59) 



While it will be convenient to use our initial definition of the coherence functions in practical 
calculations, we will always present our results using the canonical normalization. 

The fact that fC c h± can be expressed as products of covariant projectors onto opposite 
energies allows a nice pictorial representation of the propagators shown in figure 1: the 
mass-shell propagator K. m h + corresponds to the usual positive energy states moving along 
the direction of momentum k and /C m h_ to negative energy states that move against the 
direction of k (that is, their physical momentum is — k). The coherence propagators on the 
other hand, can be viewed as transporting states simultaneously in from both directions 
{tCch-) or out to both directions (JC c h+), both along and against the direction of k. From 
practical point of view however, it is important that the matrices IC m h± and lC c h± form an 
orthogonal system with respect of Dirac traces. Indeed, using first the definitions Eq. (2.53) 
the collision integrals (2.38-2.39) can be written as: 



Cmh± [fa] 
Cch± [fa] 



1 

2w k 

A. 

2w k 



Tr 
Tr 



(Ccoii + C ( | oll )7 /C m h± 
(Ccoii + C^ oll )7°/C cft=F 



(2.60) 



where £ k = ^ k /(^ k ~~ m i?)- Now using the orthogonality properties of fCjh± and the generic 
expression for the collision term Eq. (2.49) we find 



Cmh± [fa] 
Cch± [fa] 



[ ^mh±fmh± + ^chT^chT 



> <H> < 



,(2.61) 



where the functions £*^± are the following traces of the self-energy functions S^^: 



<>>. 



J mh± 



y<,> 
^ch± 



1 



2w k 
1 

2o; k 



Tr 



-Tr 



iS dt > (± w k)£mfct ( k ) 



iX</(±u; k )/C c/iT (k) 



(2.62) 



All other contributions from the collision integrals {C co \\) get annihilated by the projection 
operators inside the traces. Note that the coherence-shell self-energy functions get 
projected with an "inverted" spin structure K, c h^ under the trace. If the coherence func- 
tions are set to be identically zero, and the source terms are neglected in the flow terms, 
equations (2.35) with the collision integrals (2.61) reduce to the usual Boltzmann equation 
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for the mass-shell excitations. Observe however, that any departure from equilibrium, even 
in the absence of the source terms, would lead to creation of coherence through £k-terms. 
Coherence would be created in this way for example in massive particle decays, which could 
be of relevance for example for leptogenesis calculations. 

The quantum Boltzmann equations (2.35-2.36) together with the explicit forms of the 
collision integrals Eq. (2.61) constitute one of the main results of this paper. We wish to 
stress that they are generic equations in the cQPA framework, where arbitrary couplings 
of the coherently evolving system to its surroundings are encoded into the perturbative 
expressions of the self-energy functions (2.62). As we saw in section 2.4.1, to compute 
these functions we need special techniques that account for the resummations over the 
coherence propagators. We will develop the necessary calculational rules in section 4, after 
we first extend the results of this section to the case of scalar fields. 

To conclude this section, let us stress that despite their apparent completeness, the 
cQPA quantum Boltzmann equations (2.35-2.36) correspond to an approximation scheme. 
Indeed, we found the spectral representation for the correlator in terms of the eight f a - 
functions by expanding and solving the collisionless constraint equations (2.9) to the zeroth 
order in gradients. Because these equations were purely algebraic they had a formally ex- 
act spectral solution which was in one-to-one correspondence to the lowest moment of the 
correlator S. If we included gradient corrections into the constraint equations however, a 
more complicated set of distribution functions with new independent terms (for example 
like ~ aa^dk S(k 2 — m 2 )) would be needed to find closure consistently order by order. In 
this sense the qBE's (2.35-2.36) written in terms of / a 's are almost too promising; in reality 
the heart of the cQPA-scheme remains to be that it allows a reasonable approximation for 
the collision integral in the master equation (2.27) for the lowest moment of the correlation 
function S. Adding more gradients to constraint equations would lead to a more compli- 
cated ansatz, and correspondingly, to find a complete set of qBE's with the new associated 
shell-functions, more moments of the dynamical equations would be needed. Most phys- 
ical observables can be written in terms of the lowest moment functions however, and in 
most physical applications the gain in the accuracy from going to higher orders in gradi- 
ents would be very limited. Our restriction to the lowest moment expansion (lowest order 
in gradients) is, at any rate, exactly analogous to the derivation of the usual Boltzmann 
equation. 

3. cQPA for scalar fields 

We now formulate the cQPA formalism for scalar fields using analogous approximations to 
the ones we introduced for fermions in section 2.1. As was shown in [13], one obtains a 
qualitatively similar phase space structure, with the mass shells at ko = ±Wk and coherence 
shells at ko = 0. However, the integration procedure leading to a closed set of Boltzmann 
equations is somewhat different from the fermionic case. First, since there is only one 
component in the scalar field correlator A < (k, t), one needs to introduce a finite number of 
/co-moments of the initial singular correlators [13]. Second, as the constraint equations for 
scalars are not algebraic there are no formally exact spectral solutions to them. In other 
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words, the gradient expansion is slightly more delicate for scalars than for fermions. As a 
result, it will be more convenient to present the scalar qBE's in terms of moments and use 
the spectral solutions only to compute the collision integrals. 

We begin with the Kadanoff-Baym equations for the Wightman functions 

iA < (u,v) = (v)</>(u)) and iA > (u, v) = ((f>(u)^(v)) . (3.1) 

In the Wigner representation, the Kadanoff-Baym equation for these functions become (see 
e.g. ref. [5]): 

(k 2 - \d 2 x + ikod x - m 2 e~i A<'> - e^{n H }{A<>>} - e~ l< >{Il<>>}{A H } = C coI1 , 

(3.2) 

where A H = A F - (A> + A<)/2 and U H = U F - (IP +U<)/2, while A F and U F denote 
the time-ordered (Feynman) propagator and the corresponding self-energy. The collision 
term can be found from Eq. (2.5) with S — > A and S — > II. Proceeding through steps 1-3 
in the approximations detailed in section 2.1 and breaking equations into hermitian and 
anti-Hermitian parts we find: 

(k 2 - \d 2 - cos (\d?d^m 2 yA<>> = -C A (3.3) 
(k Q d t + sin{^d?d£ )m 2 )iA<>> = C H , (3.4) 

where C H = (C coll + C c f oll )/2 and C A = (C coll - C c t oll )/(2i). 
3.1 Phase space structure 

We now proceed to approximation step 4 and analyze the KB-equations (3.3-3.4) in the 
zeroth order in T and dtm to find out the singular phase space structure. That is, we 
initially set: 

(k 2 - ^ - m^iA^ = 
/c a^A < - > = . 

Because both of these equations contain explicit time derivatives even in zeroth order, 
one need to use them both to get one algebraic constraint. The appropriate approximate 
solution was found in [13] : 

*A<>> (k, t) = 2n sgn(fco)/£> (|k| , t)5(k 2 - m 2 ) + 2n /<•> (|k| , t)5(k ) , (3.5) 

with Sfc = sgn(/co), and m = m(t). Following the fermionic analog, it will be convenient 
to define two new (dependent) coherence-shell solutions: 

f<±>=(u k ± l -d t )f<>>. (3.6) 
With these variables we can write Eq. (3.5) as: 
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iA<>>(M) = — E (±/™± > (k,^(fcoT^) + / c f(k,i)^o)) • (3-7) 

As in the fermionic case, the KB-equation for the spectral function A = |(A > — A < ) 
is identical to the ones for A <>> , and consequently, the solution is of the same form as 
Eq. (3.7). In addition however, the spectral function must obey the sum rule, which 
follows from the equal time commutation relations of the field operators <j>: 

d ^(k + i -d t )A(k,t) = l. (3.8) 

7T 2 

Again the spectral relation completely determines the spectral on-shell functions, setting 
f±=\ and f^± = (see ref. [13]), reducing A to its standard form: 

A = ir sgn(k )5(k 2 - m 2 ) . (3.9) 

Using this result with the defining relation 2iA = A < — A > , one can easily show that the 
dynamic functions / > and f < are related: 

/> ± = l + /< ± , and fe = fe. (3-10) 

That is, only half of the on-shell functions appearing in A <!> are free variables. In what 
follows, we derive equations of motion for the on-shell functions f a = /<. Finally, let us 
write down the cQPA pole propagators which are equivalent to the standard expressions: 

r,a ^ ' ^ k 2 — \m\ 2 ± zsgn(fco)e ^ ^ 

From these one can find out the cQPA Feynman and the anti-Feynman propagators: 

A F = A r + A< 

A^ = -A a + A<, (3.12) 
where A < is of course given by Eq. (3.7). 
3.2 Equations of motion 

We again define the cQPA transport equations by treating the functions f a as free param- 
eters, and inserting the spectral solution (3.7) as an ansatz back into the full KB-equations 
(3.3) and (3.4). Because we have only one scalar function A < for three unknown shell func- 
tions, we need to integrate these equations with a number of different weights functions to 
get a closure. To be specific, we will use the moment functions: 

Pn (k,t) = J ^AffiA<(M), (3-13) 

which we need at least three to get the closure. Taking the three lowest moments and 
working to the zeroth order in gradients, we find the following invertible relations between 
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the moments and the on-shell functions f± 



and /<: 




Pi 



\{U + f-) 



(3.14) 



We will clearly need three evolution equations for our three moments, and it is natural to 
define them as the zeroth moment of equation (3.3) and as first and second moments of 



It is now evident that the moment connections (3.14) are the key element of the cQPA 
approximation, as they will allow us to rewrite the collision integrals appearing in Eq. (3.15) 
in terms of the moments Po,i,2 to get the closure. Unlike in the case of fermions, we do 
not attempt to rewrite these equations in terms of /-functions, although it could be done 
formally by differentiating the inverted equations (3.14) and using recursively the evolution 
equations (3.15). This change of variables carries a delicate issue related to the gradient 
expansion, however. Remember that the moment equations (3.15) are exact in the sense of 
gradient expansion, and to avoid introducing a loss of accuracy in going to /-variables, one 
should treat also the relations (3.14) as exact, despite the fact that they were computed 
only to the lowest order in gradients. The inverted equations would thus have a mixed set of 
gradient terms, coming from both the exact moment equations and from the approximate 
inversion process. In particular, second order gradient terms would be invoked because 
of the second time derivative in the equation of motion for pq. These terms should not 
be neglected in the spirit of gradient expansion, however, as they arise from the defining 
relations (3.14). This issue did not arise in the fermionic case where the equations of motion 
include only first order time derivatives. As a result of this complication we prefer to write 
our qBE's in terms of the moments rather than the /-functions in the scalar case. 

However, also the scalar coherence solutions are rapidly oscillating and thus the scalar 
collision terms need to be resummed with respect to these oscillations. We can see this by 
solving the evolution of the coherence-shell solutions to the lowest order in gradients from 
equations (3.15) in the context of the formulae (3.14) and the definition (3.6): 



where the correction term O has the same meaning as in the fermionic case. Indeed, 
equations (3.16) are identical to the zeroth order limit of fermionic qBE:s (2.36). The 
mass-shell solutions on the other hand, are constants to the lowest order: dtf m ± = C 1 , 



Eq. (3.4): 



-dfpo + ulpo - p 2 
dtPi 

dtP2 ~ ^d t {m 2 )p 



- (Ca) 
(Ch) 
(k C H ) ■ 



(3.15) 



d t fc± = T2iw k / C± + O 1 



(3.16) 
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and so, analogously to equations (2.41), we find the following lowest order equations for 
the singular shell-solutions: 

d t A< ± = O 1 , d t A< ± = T2*u, k A< ± + O 1 . (3.17) 

3.3 Resummed scalar collision term 

We need to express the collision integrals appearing in Eq. (3.15) in terms of the distribu- 
tion functions f± and f c ± (and eventually in terms of the moments using the connection 
Eq. (3.14)). The basic quantity we encounter is: 

{Ca) = J ^*S^- i O({n>(fc > t)}{A<(fc > t)}-{n<(fc,t)}{A>(fc,t)}). (3.18) 

where a = 0, 1. Again, to re-sum the oscillatory gradients of the distribution functions in 
the 0-expansion, we write Eq. (3.18) in the two-time representation: 

(C a ) = \J dw (- i d ro r[u>(t+ r ^,wo)A<(w ,t- r f) 

-U<(t+ r ^,w )A>(w ,t- r ^)\ . (3.19) 

Z Z J ro=0 

For a = this immediately reduces to a formula analogous to Eq. (2.43) for fermions. For 
a = 1 the extra d ro -derivative gives rise to an additional complication, but we still continue 
to search for a consistent expansion around the external time t as before. We begin by 
writing the spectral propagator (3.7) in the two-time representation: 



iA<(«*,t4k) = j ^e-^(-o-^) A < (fco)k) ^±^o 



) 



2^k ± L J 

Here we Taylor expanded A<(k, w °+ w o ) around the external time variable t, and used the 
recursive zeroth order equations of motion: 

d?A< ± = O l , d?A< ± = ( T 2^ k ) n A< ± + O l , (3.21) 

exactly as in the fermionic case. Now, using the expanded propagator (3.20) in the collision 
integrals (3.19) we get: 

(Ca) = £([± *n> Q± /< ± + iU> a± f< ± ]-[><+<]), (3.22) 

where the effective self-energies are: 

n<'> ± = (±u, k + ld t ) a n</(±a; k ) and l£> = (ift) a D<j>(^) (3.23) 

with 

II^fco.M) = d^oe^^-^n^^^o.k). (3.24) 
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Note again the simplicity of these results; to get the collision integral for an arbitrary 
moment equation, we only need to evaluate one generic function n^ > (±Wk) ) which is of 
the same form as the fermionic effective self-energy function Eq. (2.50). For the cases 
a = 0, 1 appearing in (3.26) we have: n<' > = II<'> = n<£, II<f ± = ic^II^ + i$II<£ 
and = |d t n^, so that 

(Co)=-^E([ m eff±(/c± ± /m ± )] -[>♦+<]), 

< Cl > = \ E ( [ ^n> ff± /< ± + (^*n> ± ) ( f< ± ± f< ± ) ] - [ > o < ]) . (3.25) 

The collision integrals appearing explicitly in Eqs. (3.15) are just real and complex parts 
of these expressions: 

(C H )=3t(Co), (C A ) = $S(Co) and (k C H ) = »(Ci). (3.26) 

Equations (3.23) deserve some comments. First, the resummation of leading oscillatory 
terms in the 0-expansion pushed the (zeroth moment) coherence self-energy functions onto 
mass-shells, analogously to what happened in the fermionic case. This again conforms with 
our expectation that collisions cannot be sensibly defined for pure coherence; although co- 
herence solutions do live on ko = in the phase space, their effect is, after resummation, 
only felt as a modification of the collision rates for the mass-shell states. Second, the 
resummation effectively transforms the moment function k§ in equation (3.18), into an 
operator (ko + |<9t) Q which, after projection onto mass- and coherence shells, gives rise to 
expressions (3.23). Obviously, a naive truncation of the diamond expansion would have 
missed the derivative terms in (3.23). Note that these gradients cannot be ignored due to 
the same reason that led us to perform the diamond expansions in the first place; when- 
ever the gradient acts on any coherence function occurring in the perturbative expansion 
for n^ > , the result is not formally suppressed by mass- gradients. Thus we have to use 
Equations (3.25) as such in the collision integrals (3.26). 

In practical calculations it is most convenient to express the collision integrals in terms 
of fa's and use the inverse relations of Eq. (3.14) to write them in terms of the moments 
Po,i,2,- For the record we write the necessary inverse relations explicitly here: 

fm± = Pi ± — 92 

f c ± = (oo k ± l -d t )po - —92 ■ (3.27) 

Note the appearance of the quantity dtpo, which appears as a "hidden" variable in the 
evolution equations (3.15). 

4. Momentum space Feynman rules 

We now derive generalized Feynman rules for computing the effective self-energy functions 
S c fj and n e g through perturbative techniques, including the coherence effects. Standard 
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methods, such as the 2PI formalism (see e.g. [5]), exist for diagrammatic expansion of the 
two-time self-energies T> ab (u,v) and H ab (u, v) appearing in equations (2.50) and (3.24), and 
our task is to reduce the computation of the diagrams generated by these methods into a 
set of momentum space Feynman rules. We derive these rules as usual by replacing the 
propagators in an arbitrary diagram by our resummed propagators (2.48) and (3.20) and 
performing all time integrations related to internal vertices. The only essential complication 
comes from the nontrivial phase structure associated with the coherence-shell parts in the 
dynamical propagators (2.48) and (3.20). We can account for all phase factors by rewriting 
all propagators in a generic 4-dimensional representation as follows: 

GW w' , k) = / ^ e -iko(wo-cw> )+ik (i-c)t G ^ k t) _ (41) 

J 27T 

Here G can refer either to fermion or scalar propagators (the internal degrees of freedom, 
included in S(ko,'k,t), are not necessary for our treatment here): G = S r , a , S^± , Sccft± 
or G = A rja , A^'±, A^g ± , where are defined below in Eq. (4.2). The sign fac- 

tor c = c(ko) has the "normal" value c = +1 for the dynamical mass-shell propagators 
G = Gm > and for the pole propagators G = G r>a which do not contain any coherence so- 
lutions. However, for the resummed coherence-shell propa gators G = the sign factor 
is negative, c = —1, as required by Eqs. (2.48) and (3.20). In the latter case the overall 
phase is "abnormal", except in the particular case of w' = t, where the c-terms cancel, 
and the phase factor becomes normal also for the coherence propagator. This implies that 
the phases associated with a given vertex in a self-energy diagram are normal in all cases 
but those where the vertex time (not equal to t) corresponds to the second time argument 
of at least one coherence propagator connected to the vertex. In these cases the signs of 
the phases coming from the corresponding coherence lines are reversed and extra phases 
proportional to the reference time t are added. This is a general rule to be used in addition 
to the usual combinatorics after the interactions have been specified. 

Let us stress that the effective mixed representation propagators Gr C)e ff(&o, k, t) in (4.1), 
corresponding to the resummed coherence propagators (2.48) and (3.20), by definition have 
their poles on the mass shells. Let us recall their explicit expressions: 

iS%> = iS<£ + iS<>> ± = 2n(S<£ + S%>) 5(k T c k ) 

iA<£ = zA<<> + *A c <;> ± = ± f<£ + /<;>) 5(k T ^ k ) , (4.2) 

where S- matrix functions are defined in Eq. (2.26). Let us remind that in terms of the 
complex path indexing the Wightman functions correspond to the off-diagonal propagators 
S* 12 = S K and S 21 = S* > . In loop calculations we often encounter also the diagonal Feynman 
and anti Feynman propagators S" 11 and S 22 . Because the novel coherence solutions only 
appear in the Wightman functions, we can find the coherent S^-functions directly by using 
Eqs. (2.24) and (3.12) together with Eqs. (4.2) and the expressions for the standard pole 
propagators given in Eqs. (2.22) and (3.11): 

oil c c< c 1 C~ 
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Figure 2: Generic diagrams contributing to the fermionic self-energy function E e g (left) and to 
the scalar self-energy function IT c ff (right). 
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where Sf,o, &f,o, S f and refer to the standard vacuum Feynman- and anti Feynman 
propagators. Moreover, the quantities S~ s and A^ correspond to the effective Wightman 
functions and A< s from which the vacuum parts have been subtracted off. Indeed, one 
can show for example that 

S r = S Ffi + 2mQ{-k )($ + m R - i^m^k 2 - \m\ 2 ) , (4.4) 

and when the delta-function is absorbed into S^ s it exactly cancels the vacuum part (in- 
cluded in the negative frequency mass-shell function S*<_ oc f<_ = 1 — n) from S~ s : 

S cK = S cs ~ 2TriO(-ko)Qk + m R - i-/ 5 mj)5(k 2 - \m\ 2 ) . (4.5) 

It is straightforward to show that similar relations hold for all propagators in (4.3). 

This completes our rules for the propagator functions in the 4-dimensional mixed state 
representation. Before we turn to the derivation of the vertex rule, let us write down the 
Hermiticity properties of our propagators in the mixed representations 

(iG{k,tyf = iG(k t t) 
(iG<j> ± (k,0 =iG<j>^(k,t), (4.6) 

where iG = , S r>a , iA m ' or iA. rA in the first line and iG = S or iA in the second. 
Note that taking the complex conjugate flips the ± — )• =p in the coherence propagator. 
Finally, it will be useful to observe that for a real scalar field propagator 



Ar,a(fc)t) — A 0)r ( k,t) , 
A> i± (M) = A< i± (-/c,t), 

A> e ff±(M) = A< eff± (M) = &c,eS±(k,t). (4.7) 

Note in particular that there actually exists only one distinct type of coherence propagators. 
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Figure 3: A cut through an arbitrary diagram around the continuous fermion line contributing to 
S c ff. The scalar lines may be interconnected or connected to separated closed fermion loops not 
shown. 



4.1 Vertex rules 

In order to define vertex rules we need to specify a concrete model. We shall adopt the 
following Yukawa interactions between our scalar- and fermion fields: 

Ant = -y^^tp + h.c. (4.8) 

Generic examples of self-energy functions generated by this interaction are shown in fig- 
ure 2. Let us consider the fermionic self-energy diagrams first. The most general diagram 
has one continuous fermion line connecting the initial and final times, accompanied by an 
arbitrary number of closed fermion loops where all fermionic lines can be connected and 
associated with an arbitrary number of scalar lines. We wish to derive a generic mixed 
space representation for such a diagram and extract the local Feynman rules from the re- 
sulting expression. In particular we will need to show that the global phase proportional 
to the external time t, arising from our propagators (4.1) vanishes. 

Consider first a generic sub-diagram of the type shown in Fig. 3. This is the continuous 
fermion line going through an arbitrary fermionic self-energy graph. It contains n vertices 
and n — 1 fermion propagators, and the final time on the fermionic line is the special 
external time t. The scalar lines may either be interconnected, or they may be connected 
to closed fermion loops, not shown in the diagram. To be able to associate correct phases 
with the vertices, we first have to introduce the notion of the ordering of the arguments 
in the two-time propagator (4.1) into the mixed representation. This can be done by 
associating propagators with a unique direction of flow. We start by defining the flow 
in the two-time representation according to the flow of time in time-ordered propagators 
G 11 (it, v); that is, from the vertex v to the vertex u. This choice induces a natural definition 
for the flow direction in the mixed representation as being along the four-momentum of 
the positive energy states. For fermions this corresponds choosing the direction along the 
fermion number flow. Given this definition, Eq. (4.1) implies that a fermion propagator 
with an energy g^o, connected to a vertex Wjo gives rise to a phase factor that depends on 
the direction of the flow according to: 

e -iq ko w j0 "incoming" 

e ic q kqkoWjo+iqko{i-c q k)t "outgoing" . (4.9) 

For a neutral scalar line we do not have a similar natural orientation 5 and we have to use 
This is so because we are considering a neutral scalar field here. For a charged scalar field the charge 
conservation could be used to define a natural flow orientation. 
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the general expression: 

e isjC vk p ko Wjo+ip ko \(l+Sj){l-c pk )t ^ (4-10) 

where the index Sj = 1 for an outgoing and Sj = — 1 for an incoming scalar line. The 
definition (4.10) is clearly consistent with the fermion phases in Eq. (4.9). (Note the 
implicit rule that for Sj = — 1 we always have c p k = 1 in the associated propagator.) 

With these definitions it is easy to show that the time integration in each of the internal 
vertices in the diagram in Fig. 3 (j runs from 1 to n — 2) gives rise to a delta-function: 

J dw 0j e^jo^o+^wjo^o+ijo)^* = 2vr S(c qj q j0 + s jCpj p j0 - <Z(,+i)o)e l<M , (4.11) 

where the extra phase factor is 

<f>j = (1 - c qj )q j0 + * (1 + sj)(l - c pj ) Pj0 . (4.12) 

Clearly the extra phase vanishes if c q j = 1 (normal fermion propagator) and either Sj = —1 
(incoming scalar field) or Sj = 1 and c p j = 1 (normal outgoing scalar propagator). Taking 
into account the special vertices at the ends of the fermion line, and integrating over the 
times at internal vertices as well as over the specific time wq, the generic expression for the 
self-energy contribution from diagram 3 becomes: 

£</ (fc ,k,i) = J dw e lk ^- W0 ^(t,w ,k) 

n 

x ( II 27r6 ( c QjQjo + s J c PjPjo ~ %+i)o))e^ linc * x C* loops , (4.13) 
i=i 

with the understanding that q n o = ko, and the global phase factor is 

^ n— 1 

<Aine = k - q w + s c p0 p 00 + - (1 + s )(l - c p0 )p 00 + ^24>j- (4.14) 

i=i 

In the second line of Eq. (4.13) we suppressed all momentum integrations and the effective 
propagators associated with the internal lines. They can easily be inserted back afterwards. 
The explicitly shown delta functions and phases come from the vertices in the diagram 3 
and the factor Ci oops contains the contributions from all possible closed internal loops in 
the full diagram. By a recursive use of the delta-functions in the internal vertices, one can 
show that (f)\- mc can be associated with a sum of scalar momenta: 

<Aine = ^(1 + S )P00 ~ ~(1 - S )Cp0P00 



n—1 



+ k - q 20 + siCpipio + -(1 + si)(l - Cpijpw + ^2 

n-1 



J'=l 

n-1 
3=0 
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Figure 4: A generic closed fermion loop diagram with all internal vertices. The scalar lines can 
again be interconnected, or connected to some other loops or to a continuous line going through a 
fermion self-energy diagram. 

where in the last step we used the fact that for an incoming particle with Sj = —1 we 
always have c p j = 1, and that for an outgoing state with Sj = 1 the c P j-term vanishes in 
Eq. (4.15). The final sum in Eq. (4.15) would vanish if all scalar lines in the graph were 
interconnected (no closed fermion loops in the graph), because then each energy piQ would 
appear twice in the sum, once with Sj = +1 and once with Sj = — 1. However, as some of 
the lines may be connected to loops, the sum does not vanish in general. Physically this 
means that coherence information can be transported between separate fermion lines by 
the scalar fields. 

Now consider a closed internal fermion loop depicted in Fig. 4. For a closed loop the 
direction of flow could be chosen arbitrarily, but we follow our earlier definition of going 
along the fermion number flow. The calculation proceeds similarly to the case with the 
open fermion line, the sole difference being that now all the propagators and vertices are 
internal ones. After a straightforward calculation one finds that for a loop with n% vertices: 

m 

Cioop « II 27r< H c i9jo + SjCpjPjo - g (i+1)0 )e l ^W (4.16) 
i=i 

where 

^loop = Yl $3 
3=1 

m 

= qil - Q20 + + Sl)piO - |(1 - Sl)Cplp j0 ] + ^ (pj 

i=2 

rii 

= QlO - C qn q n - S n Cp n p n o + ^ [|(1 + Sj)pj - 1(1 - Sj)c p jPjo] 

3=1 

rij 

= Y2s jPj o, (4.17) 

3=1 

where in the last step we used the fact that due to cyclicity n + 1 'th and the first fermion 
propagators are the same: c nq q n o + s n c pn p n o = (?(n+i)o = QlO- Combining the phase factors 
from the open fermion line and all closed fermion loops we find that the total extra phase 
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Figure 5: Generic diagrams contributing to the scalar self-energy function II G ff where the special 
times t and wq are connected with the same (left) and with different (right) closed fcrmion loops. 

proportional to the external time t, counting all 2N vertices in a diagram with a total of 
N internal scalar lines is: 

27V 

<pTOT = 4>line + £ $oop = £ S lPlO = • (4.18) 
i 1=1 

The total phase </>tot vanishes because each scalar propagator appears twice in the sum, 
both as an outgoing (s/ = +1) and an incoming (si = —1) one with the same energies, 
and these contributions cancel pairwise. Physically the vanishing of </>tot m Eq. (4.13) 
shows that while the energy is not conserved in the internal vertices in connection with the 
coherence propagators, the overall energy in the complete diagram is conserved. 

Showing the vanishing of the total phase factor for an arbitrary scalar self-energy 
function Il e ff proceeds similarly to the fermionic case. In addition to the closed internal 
fermion loop one now has to evaluate the phases coming from the two sub-diagrams shown 
in Fig. 5, which express the two possible ways of connecting the special vertices t and wo 
into a most general diagram for II e ff (t, wq, k). It is by now quite straightforward to show 
that the extra phase, in either of these cases becomes 

n— 1 n— 1 

<Afigs.5 = ^2 fa + ( X ~ C 9«)*0 + <7n0 - <?10 + k = ... = ^SjPjQ . (4.19) 

Combining this result with the all possible phases from the internal closed loops gives the 
total phase factor for a II e ff with N internal scalar lines: 

27V 

</>TOT = 0figs.5 + £ 4op = E = > ( 42 °) 

i 1=1 

where the final sum over the scalar energies vanishes by the same argument as in the 
fermionic case in Eq. (4.18). This result completes our proof that the local extra phases 
coming from the coherence propagators cancel in arbitrary self-energy diagrams. As a 
result, we can neglect all such phases in the actual calculations and use the local momentum 
space vertex Feynman rule: 



(fyiptf) : y {2ti) a 8{% - c g q - s p c p p )S 3 (<i - q - s p p) 



(4.21) 



-ij/(27r) 4 5 3 (q'-q + p) 

xS(q' -c q q +p ) iS$(q,t) 



b 



-ly^YS^ - q - p) 

X<%£ - C q q ~ Cj,po) « A eff(P^) 
a P b 



Figure 6: The Feynman rules including coherence with an explicit orientation of the scalar lines 
following from Eqs. (4.2 - 4.3) and Eq. (4.21). The sign functions are c P:q = —1 for the coherence 
parts of the associated propagators and c p . q = 1 otherwise. 

The only difference from the usual rule then is the appearance of extra sign factors in 
the energy delta function in association with outgoing coherence propagators, leading to 
a local energy non-conservation within the loop. The complete set of momentum space 
Feynman rules for computing the fermionic and scalar self-energy functions including co- 
herence propagators in the Yukawa theory with interaction Eq. (4.8) are shown in figure 
6. The arrows in the propagator lines indicate the direction of the flow corresponding to 
the 4- momentum of the positive energy state. 

5. Examples 

We shall now compute some examples of self-energies and collision integrals in the Yukawa 
theory described by the Lagrangian (4.8), starting with with the simplest one- loop self- 
energy diagrams shown in the upper right panel of Fig. 7. These diagrams can be obtained 
for example from the two-particle irreducible (2PI) effective action vacuum diagram shown 
on the left in Fig. 7: 

r 2 pi = --£/" d\d 4 vTr{S(u,v)S(v,u)} A(u,v) , (5.1) 
2 Jc 

where the integration is over the Keldysh path [11]. For example the fermion self-energy 
now follows by a direct functional differentiation: 

£ afe (u, v) = = S ab (n, v) [A ab (u, v) + A ba (v, u)] . (5.2) 

This much of the calculation is straightforward even without the Feynman rules developed 
in the last section. However, to proceed further would be tedious, since the two-time 
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Figure 7: On left: the diagram contributing to the 2PI effective action at the lowest order for 
interaction (4.8). On top right: the two diagrams contributing to the fermion self-energy £"fj(fc,i) 
and on bottom right: the single diagram contributing to the scalar self-energy II^(fc,i). 



representation of Eq. (5.2) should be integrated over wq according to Eq. (2.50) using the 
two-time effective propagators given by Eqs. (2.48) and (3.20), while taking great care of 
the different phase factors in the number of coherence- and mass-shell propagators. 

With the Feynman rules of Fig. 6 at hand none of the extensive labour discussed 
above is needed and we can directly write down the final mixed representation expression 
for T, eS (k,t): 

y 2 sr^ f d 4( 7 d 4 p 



i^{k,t) = V - J (2vr)4 ( t*4 «ffig'*)P A gfo*)*'.i + <(P.*)Vil >< 

Sfc — i 1 

x (27r) 4 5(/c - c q q - c p [s b ]po)5 3 (k - q - sp) . (5.3) 



Three observations are in place here: first we used a sum over the index s b together with the 
Kronecker delta functions associated with the scalar propagators to account automatically 
for the direction of the four momentum flow. Second, we defined G® b to denote effective 
propagators Eqs. (4.2 - 4.3) for which the coherence parts G^f'g are combined with the 
correct c Pi fc-factors in the vertex 5- functions. Finally we took care of the implicit s b - 
dependence of c p -factor in the scalar line by introducing the notation 

c p [s b ] = c p 5 s b :l - 5 s b _! (5.4) 

inside the vertex delta-function. These notations will be useful later on. Note that the 
implicit dependence in c signing out the coherent parts of the propagators still remains as 
explained above. 

Expanding the sum in s and separating the coherence solutions, we find that for ex- 
ample the self-energy functions iT l < = iS 12 and iT, > = iS 21 become: 

z£</(M) = 2/ 2 f ( ^)4 ( ^ 4 (2vr) 3 ^ 3 (k- q -p)x 

(2vr)(5(A:o - <?o - Po)iS% > (9^)«A<' > (p, t) 
+ (2Tr)6(k + qo-poliS^fafy^fat) 
+ (27r)<5(fc - qo+p )iS<' > (q,t)A CfiS (p,t) 
+ (2ir)5(k + q +p )iS<>>( y q,t)A c , cS (p,t)] , (5.5) 
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where we used the fact A> (— p, t) = A< (p, t) and A Cje ff (po, — p, t) = A C)C s(pq, p, t) (because 
of the isotropy and the identity A> eff = A^ cff = A CjC g). Equation (5.5) shows explicitly 
how different types of energy conservation are associated with the coherence propagators; 
the delta-function associated with the 5 m A m -term has the normal signature. Using the 
isotropy again we can combine the coherence and mass terms together under the same 
delta function 5 A (k — q — p) which can then be used to integrate over the momentum p, 
eventually giving just: 

l ^(k,t)=y 2 J ^4</( 9 ,i)iA</(fc -<?,*), (5.6) 
where we have further defined 

iA<£ = ^(±f<!> + f<>) 5(k T u*) . (5.7) 

Note that these functions differ from the effective propagators in Eq. (4.2) in that here the 
coherence-shell functions appear in "wrong" energy shells. 

Equation (5.6) is remarkably simple; it can be obtained from the standard expression 
for S <>> by a direct substitution S < ' > — > S^f* and A < > > — > A^ > . Unfortunately such a 
simple rule does not generalize to arbitrary diagrams, as can be seen already from the one- 
loop scalar self-energy function n e ff • Indeed, a direct evaluation of the scalar self-energy 
diagram shown in Fig. 7 gives 

•n&(M) = -y / ff)^ t^ [iS t ^ iS 4 fo> *)1 x 

x (27r) 4 % 20 - c\q w - fe )<5 3 (q 20 - qio " k) , (5.8) 
from which it is now easy to see that the self-energy functions n^ > become: 

tn%>(k,t) = -y / j^TtiiS^fatyS^ik - q,t)] • (5.9) 

That is, the substitution of S < ' > — > S^ > is made only for the fermion line which is flowing 
out from the special vertex t, while for the fermion propagator flowing into the vertex 
the substitution is S <,> — > . These simplifications generalize for propagators directly 
connected to the special t- vertex in arbitrary diagrams. In all other cases the assignment 
of energy signs in the coherence propagators is more complicated and can only be worked 
out by a use of the full Feynman rules of section 4. 

5.1 A two-loop example 

Let us next consider a more complicated 2-loop example. The contribution from the 
diagram shown in Fig. 8 to T> ab is: 

fe ,e ,/ ,6 J i=l K ' j=l V ; 



- 30 - 




Figure 8: The two- loop 2PI-fcrmion self-energy diagram. Each scalar propagator goes in both 
directions, so the diagram represents four independent terms. The special vertex a corresponding 
to the external time t is again marked by a black dot, while the internal vertices are shown by open 
dots. Full explanation in text. 

X (2vr) 12 5 3 (q 2 - qi - 4 P 2)<5 3 (q 3 - q 2 - ^ Pl )5 3 (k - q 3 - s b p 2 ) 

X 5(q 2 ~ C ql q 10 ~ C^ 2 p 2 o)%30 - C 9 2<720 - CpiPlO^Oo - C q mO - C b p2 P20) 

x [iA^( Pl )5 shl + iAS(pi)<y a ;,-i] [iA b J p2 ( P2 )5 s f A + iA{ b 2 {p 2 )5 sb J , (5.10) 

where we have suppressed the i-arguments in propagators for clarity. We also continued 
using the tagging on propagators G ab h introduced in the 1-loop example above, as well as 
the Kronecker delta notation and s^-sums to indicate all possible directings of the scalar 
vertices. We also used the shorthand dp = c p [Sp] for the notation Eq. (5.4) in the vertex 
delta functions. These notations unambiguously indicate the scalar momentum flow in the 
internal vertices and whether the coherence solution is to get a nontrivial sign or not. 

Now, because the momentum pio appears only in one of the delta-functions, we can 
absorb the index c e pl to a scalar propagator by effecting a transformation pio — > c p i[s e ]pio, 
which amounts to the substitution A^(pi) — > A"|(pi). The same argument applies to 
the energy variable qio corresponding to the fermion propagator flowing into the special 
i-vertex, which we can deal with by a change q\o — > c q iqio, causing S^i (<7i) — > S^(qi). 
Beyond these simplifications associated with the special vertex t no significant reduction is 
possible. Being careful in particular with the scalar flow direction assignments on vertices 
e and / one eventually finds: 

x zS^q^Sfjq^Jjp,) + iA{ b p2 (-Cp 2 p 2 )) 

x (27r) 12 <5 3 (q 2 - qi - p 2 )5 3 (q 3 - Q2 - Pl)£ 3 (k - q 3 + p 2 ) , 

X %20 - 910 - C p2 P20)%30 ~ C q2 q 20 ~ Plo)^(^0 ~ Cq3<730 + P2o) ■ (5.11) 

The effective Wightman functions and A^ > were defined in Eqs. (5.7). The special 
diagonal propagators G^ and G^g {G again denotes S or A) are given by the analogous 
relations to Eq. (4.3) with Gf s replaced by Gf s everywhere, and the functions G~ s again 
correspond to functions from which the vacuum parts have been subtracted out. 
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Particular self-energy functions, and particular sub-contributions to these self-energies 
are obtained by assigning special values for the time-contour indices a, b, e and / in 
Eq. (5.11). Let us consider one example here for illustration. Choosing aefb = 1112 we 
get a two-loop correction to the self-energy function S< = S 12 : 

•E£(M)ni2 = yJ ^4 ^ ^ in (g3,P2,t)^ 2 3(g3)[^; 2 ( P2 ) H-iA^-oaa)] 

x(2tt) 4 (5(A;o - 0,3-730 - c p2 p 20 )5 3 (k - 03 - P2) , (5.12) 
where the function A 111 is given by 

zA^(q 3 ,p 2 ,t) ee y 2 | ^^L[^ fe+ p 2 ),A^( m )]^ c 1 g 1 2 fe) 

x(2vr) 4 5(g 30 - c g2 920 - Pio)£ 3 (q3 - 02 - Pi) ■ (5-13) 

Note that the self-energy function (5.12) becomes equivalent with the 1-loop result shown 
in Eq. (5.3) when iA 111 — > 1. In this sense A 111 can be interpreted as a one-loop vertex 
function correction to the 1-loop self-energy. On practical side, observe that the coherence 
parts of S^ 3 (qs) in Eq. (5.12) and S^ 2 (q 2 ) in Eq. (5.13) are associated with a different phase 
space integrals than the non-coherent parts, due to the Cgj-dependence of the remaining 
energy delta functions. This dependence is passed onto the arguments of the other propaga- 
tors when one performs the integral over the remaining delta function. Of course there are 
other contributions to S < and also S > is needed to compute the corresponding collision in- 
tegrals in the quantum Boltzmann equations. Nevertheless, we believe that these examples 
amply display the feasibility of the use of our Feynman rules for practical diagrammatic 
calculations in coherent perturbation expansions. In particular, a two-loop calculation 
similar to the one presented above, but generalized to the multiflavour case [16], will be 
relevant in the context of resonant leptogenesis [8]. 

5.2 Direction independent Feynman rules 

The result (5.10) for the first time displays fully internal scalar and fermion propagators. 
In particular, it shows how the internal scalar propagators combine to a form which is 
insensitive to the choice of flow directions of the scalar lines. We can use these findings to 
reformulate the direction dependent Feynman rules of Fig. 6 in an even more useful form. 
The complete rules can now be stated as follows: 

• Draw the diagram just as in the vacuum field theory and associate the usual vacuum 
symmetry factor with it. 

• Give each vertex an index Sp associated with the direction of the scalar field, and use 
the vertex Feynman rule 

(j)^: -iy (2ir) 4 5{q' -c q q -c p [s b p }po)5 3 (c{ -n- s b p p), (5.14) 

where c p [sp] = c p 5 s b i — S s b_i, where c p = —1 for the coherence parts of the propa- 
gators, and c p = 1 otherwise. 
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Figure 9: The final direction independent rules including the coherence. 

• For each fermion line in the diagram, substitute a sign-tagged propagator S^(q,t), 
where the index c q explicitly keeps track of the coherence sign in the effective fermion 
propagator S^(q,t) defined in Eqs. (4.2-4.3). 

• For each scalar line substitute a propagator 

S s a_ s b , (5.15) 

where A^(p, t) is again the sign-tagged propagator corresponding to the effective 
scalar propag ator A^(p,t) defined in Eqs. (4.2-4.3). 

• Sum over all indices Sp associated with all vertices. 

• All other Feynman rules associated with the momentum integrations and negative 
signs associated with closed fermion loops are as usual. 

Special simplifications of the rules apply for propagators connected to the external vertex 
t; see the discussion below Eq. (5.9). These direction independent rules are presented 
graphically in the figure 9. 

5.3 Explicit 1-loop collision integrals for fermions 

We now compute the final expressions for the fermionic collision integrals (2.61) following 
from our above examples for S e g. To keep our expressions simple we only consider the 
1-loop diagram, and that the scalar field is in thermal equilibrium with (see eg. [5]): 

iA<^(p) = 2n S gn( Po )f< q >l {po) (p)8(p 2 - mj) , (5.16) 

where /<|(p) = l/(e±^p - 1) and />£(p) = 1 + /<£(p) with f3 = 1/T . Using the 
equilibrium propagator (5.16) in Eq.(5.6), going back to the form with both q- and p- 
integrals and integrating over go and po using the on-shell delta functions in the propagators, 

-33- 




iD 



ah 



iA?(p,t)5 s , 1+ iA b Vp,t)5. 



we find 



x *( ±Wk - SqWq - apW j)/<;>f (pJEg^q) + 5 S- Sq (q)l -(5-17) 

The sign- factors s qiP appear in Eq. (5.17) as a result of projections on positive and negative 
energy shells after integrations over go and po; they should not be confused with the 
orientation signs encountered in our Feynman rules. According to the standard 3-particle 
kinematics the momentum delta functions have roots only if > 2\m(t)\ and only for 
the signatures, 5(a;k + oj q — ojp) and <5(— uj^ — oj q + uip), so that each on-shell self-energy 
function has only one contribution coming from the sum in Eq. (5.17). Employing the 
decompositions (2.53) we finally get: 

fflSf»(±<*,k.i) = 4^)3 ^fe (*>'«("■■ + - «> 3 ^ - 1 - p) 

x/S*(p)E(Wq)& > T h') + t*±(q)^(<tt)) (5.i8) 

ft' 

Given Eq. (5.18) it is a simple matter to calculate the scalar self-energy functions E^ ft ± 
which appear in Eqs. (2.61) 6 . Multiplying (5.18) by the projectors K.jh± and taking the 
traces we find: 

^ - / wife 4 '- + - "> 3{k - q - p) /snp) 

X VJ (MV ± (tf, fc, q , k)/^(q; f) + X| ± (ft', ft, q, k)/<;J(q; ()) , (5.19) 

h> 

where the generalized squared matrix elements are given by the traces: 
M 2 mj± (h',h,ci,k) = y 2 Tr[/C m ^ T (q)/C^±(k)] 

M 2 cj± {h', h, q, k) = y 2 TV[/C cV ±(q)/C j7l ±(k)] , (5.20) 

with j = m, c. Note that these matrix elements are not in general real (see Eqs. (5.24-5.25) 
below). Given the self-energy functions we can now combine all terms appearing in 

6 Note that the signature of the uiq and q in the energy- and momentum conservation delta functions are 
different. This is because one of the fermions labelled by k and q must be an antiparticle corresponding 
to a negative energy state in our language. That is if, say, the k-state is a positive energy particle, then 
the q state must have a negative frequency. However, since in Eq. (5.18) we wrote the energies in terms of 
positive physical frequencies, which makes the momentum q to appear with a "wrong" sign. The solution 
is of course that here the physical 3-momentum of the antiparticle is just — q. The same argument applies 
to the case where the k-state is an antiparticle. We could reinstate the normal signatures by a change 
k — > ±k. p — > ±p and q — > =pq, and the corresponding Feynman-Stiickelberg reinterpretation of the 
/-factors. However, there is no practical advantage of making this change and it is easier to keep using the 
unphysical momenta instead. 
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the r.h.s. of the Eqs. (2.61) to form the complete collision integrals Cjh±- The mass-shell 
collision integrals are in this way found to be: 



C m h±\fa] = J? — V / d3q N , f P (2ir) 4 5\k-q-p) X 
mh±u«\ 2^^fJ 2w 9 (27r)3 2 ^(2tt)^ ^ V ^ 

X^ m± / c < q t(p)/<,, T (q)/>, ± (k) 

+ ML±/ii(p)/^±(q)/w» ± (k) 
+ ML T /S(p)4±(q)4(k) 

+ ^^/eqt (p)/^ T (q)/ c ^(k) ] - [ >o< ] , (5.21) 

where we defined a shorthand notation 5 4 (k — q — p) = <5(wk + w q — Wp)<5 3 (k — q — p). As 
usual, we see that the distribution products in the Boltzmann equations are independent 
of the form of the interactions, the details of which are entirely encoded in the matrix 
elements. The first line in equation (5.21), which contains only the mass-shell distribution 
functions, corresponds to the usual collision integral in the standard Boltzmann equations. 
Indeed, rewriting the distribution functions in terms of the usual particle and antiparticle 
numbers, one finds for example the familiar form: 

/£t(p)/£- ("<!)/;*(*) = (1 + <(P))M<ZK(*)- (5-22) 

One could similarly rewrite all the other products of distributions appearing in Eq. (5.21) 
using the particle and antiparticle distribution functions. However, since / a 's are the 
quantities which naturally appear in all our loop calculations and in our equations of 
motion, we will stick to this notation in what follows. The second line in Eq. (5.21) comes 
from an internal coherence line in our one-loop diagram, corresponding to a coherence 
modulation of a scattering rate off a state with momentum q. The third line represents an 
on-shell contribution to the collisional coupling between the mass and coherence shells in 
the mass-shell equation, and the last term gives the coherence modulation of this collisional 
coupling term. The coherence collision integral C c h± is somewhat more complicated: 

C ch ±[fa] = — Y ! f P (2ir) 4 6 4 (k-q-p) x 

2^ k ^7 2 Wq (2^)3 2^(27r)3 ' ' 

a2 t<4> 



[ 



A44 m± / c ^(p)/<,, T (q)/> ± (k) 

+ M 2 cm± /££(pte ± (q)/£±(k) 
+ &M 2 m * c± / c < q i(p)/^ T (q)/> /i± (k) 

+^mi: ± / c < q i( P )/^ T (q)/„ > t , ± (k) 

+ (± O T)*] - [ >^< ] • (5.23) 

The interpretation of the various terms in this equation should be obvious now. Finally, 
for completeness, we give the explicit expressions of the matrix elements M. a b± following 
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from our 1-loop example. A straightforward evaluation of the traces in Eq. (5.20) gives: 



M 2 _u 

J l mm± 


= M 2 


M 2 cm± 


± m R 






M 2 mc± 


_ T m R 




M 2 CC± 


= (!-- 



M 2 ±M 2 R 



iMj 



M 2 ±M 2 R -iMj 



mR )Ml + { mR 



-)M 2 R ±i{- 



+ 



-)M 



i j 



(5.24) 



where we have defined: 



M 2 = y 2 (1 + h% q ■ k) (q- k-m 2 R + mj) + h' h\q\\k\(l - (ci-kf 



M 2 R = y 2 (1 + h'h q • k) m i? (w q - w k ) 



Ml 



(l+^q-k) m/(/i'|q| - h\k\ 



(5.25) 



where q • k = cos 6 is the angle between fermion momenta and the product of fermion four 
momenta is q ■ k = ^m^ — |m(t)| 2 . Note in particular that the squared matrix elements 
A4^ m _|_ are real and equivalent to the standard expression for a scalar field decaying to 
fermions after one makes the usual Feynman-Stiickelberg interpretation of the negative 
energy states as the positive energy antiparticles. 



6. Applications 

In this section we consider numerical examples in the case where the fermion mass changes 
abruptly but continuously from zero to a finite value. To be specific, we model the change 
by a kink profile: 

\m(t)\ = "^(l + tanh*"^ 



2 v r 
0(t) = ^(l-tanh^- t °), (6.1) 



and we assume that the fermion is interacting with a scalar 4> and another fermion field q, 
through a non-diagonal Yukawa-interaction term 

Ant = ~y QL + h.c. . (6.2) 

For simplicity we also assume that the scalar <f> is thermal and fermion q is in kinetic 
equilibrium with an effective chemical potential fi q (t). With these assumptions the self 
energies E e g reduce to the usual thermal expressions, where S <J> are related by the Kubo- 
Martin-Schwinger relation: = e ^^ 0_ ^)S < , and 

£<(£;) = (£<7°-£<k. 7^, (6.3) 
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Figure 10: Shown are the various self-energy components ^f h+ (for positive energies only) defined 
in Eqs. (6.8) as a function of trrioo for a fixed |k|/rnoo = 1 (left panel) and as a function of |k|/TOoo 
for a fixed tni^ = (right panel). For other parameters we used tjq = 0, r = 2/moo, A8 = —1.0, 



T = 10m o 



10m 



Smoo, fi q = and y = 0.3. 



where k = k/|k| and iS^g are real-valued functions, whose expressions at the mass shells 
are (remember that Wk = oj^t) depends on time through the mass function) 



iT,^(±u k , |k| 
v k iY,^(±u) k , |k| 



y 2 T" 

8vr|k 

y 2 T 2 



8vr|k 

where we defined v k = |k|/u;k and 



/i(±w k ,|k|)| 
ii(±a; k , |k| 



|J (±u; k , |k| 



/„(fe ,|k|) = 0(A) r +5 dyy n 

Ja-8 



(e»-/*«/T + l)(e*o/r-i/_i) 



(6.4) 
(6.5) 

(6.6) 



with 



|m| 2 + TTig 



2|m| 2 T 
A^dml 2 ^ 2 m 2) , k | 



2|m| 2 



(6.7) 



where A(a, 6, c) = (a + b — c) 2 — 46c is the usual kinematic phase space function. Note that 
both So and S3 have finite limit when |k| — > despite the apparent singularity, because the 
integration region oc 5 oc |k|. The chemical potential [i q is calculated from the conservation 
of the fermionic charge in the specified interaction. For the projected self-energy functions 
S^K± in the collision integrals (2.61) we find the expressions: 

1 



J mh± 



v <,> 



(1 ± / l 7J k )(iS^ > (± Wk ) " h l ^' > (±UJ k )) 



1 

2w k 



(hm R v k =F im/)(iSo ^(i^k) - hiEf^ (±co k )) 
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Figure 11: Shown is the excess total particle number density (n^ + n^) — (nk oq + rikcq) (left) and 
the magnitude of total particle - antiparticle coherence density / c (k) (right). 

Let us note that we used the self-energy function (6.3) also in the numerical examples 
in our earlier work [12]. However, at that time the relevance of the resummation of the co- 
herence oscillations was not realized, and so the coherence-shell self-energy was erroneously 
computed in ref. [12] by naively projecting T, < ' > to the ko = shell. This shows that the 
effects of resummation are nontrivial even when the primary self-energy is thermal. Typi- 
cal behaviour of the self-energies is displayed in figure 10. On the left panel we show the 
time dependence of the self-energies over the wall region, and on the right panel their 
| k| -dependence at a fixed time tru^ = 1 during the transition. In all figures in this section 
the time is measured in units and the momenta in units m^, where m,oo is the abso- 
lute value of the mass of the fermion at the infinite future moo = \m(t = oo)|. Note that 
at very early times, where the -0-field is massless, only , is nonzero as expected for 
a right chiral interaction term in association with a positive helicity state. Once the field 
becomes massive all self-energies evolve differently as a function of the increasing mass. 
In particular the imaginary parts of are very small and nonzero only inside the wall 

region, confined to |tmoo| ^ 2. In the right panel we again see that the right helicity mass- 
shell self-energy is by far the largest, while the other components rapidly become vanishing 
when the momentum increases and one again effectively approaches the massless limit. It 
should also be observed that while the mass-shell functions are strictly positive definite, 
^^+' s which appear as cross-interaction terms between the mass- and coherence-shells in 
the collision terms (2.61) can have either sign. 

We have solved numerically the quantum Boltzmann equations (2.35-2.36) with col- 
lision integrals (2.61) where the projected functions Hjh± ar e given by Eq. (6.8) for the 
mass profile defined in Eq. (6.1). Our results are shown in Figs. 11 and 12. In the left 
panel of Fig. 11 we show the evolution of the |k| -dependent excess total particle number 
density (rik + n\c) — (^keq + ^keq)> where = Y2h n kh (similarly for n), and n^h and 
nkh were defined through the Feynman-Stiickelberg interpretation in Eq. (2.17). In the 
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Figure 12: Shown is the total right chiral number density jj^k) (left) and the total charge density 
(net particle - antiparticle asymmetry) j°(k) (right). We used the same parameters as in Fig. 11. 



right panel of Fig. 11 we show the total particle-antiparticle coherence density defined as 
/ c (k) = ^ Ylh± \fch±0 Si )\, where the individual coherence functions j c h± are canonically 
normalized according to Eq. (2.59). The |k|-dependent baseline for the colour coding in 
figures was defined by the adiabatic thermal Fermi-Dirac distribution n^heq = "-kfteq = 
l/(exp(/3a;k) + 1), which ranges from 2 at the low end to ~ 1 at the high end of the plot- 
ted phase space area. Finally, we used parameters to = 0, t = 2/rrioo and A9 = —1.0 
for the kink, and T = lOm^, m q = IOtoqo and = hm^ for the self-energy. For the 
Yukawa coupling we used y = 0.3. As a boundary condition, all fields including the co- 
herently evolving fermion were are assumed to be in thermal equilibrium without chemical 
potentials in the distant past. 

Both quantities shown in Figs. 11 vanish in the equilibrium conditions, when no coher- 
ence evolution takes place. However, creation of extra excitations and coherence is clearly 
seen in the plots. As expected, most of the excitations are created in non relativistic modes 
near the transition time. As time goes on, the asymmetry diffuses towards the higher mo- 
menta in the phase space. The exact quantitative details are of course dependent on the 
kink and the interaction parameters. The creation of the extra modes seen in the left 
panel, is associated with narrow bands of damped oscillations in the coherence distribu- 
tion f c shown in the right panel. Just like the difference from equilibrium, the coherence 
function is restricted in time to near the phase transition time and in momentum to non- 
or near relativistic region. The canonically normalized coherence functions are roughly by 
a factor of five smaller than the deviations from equilibrium in the mass-shell functions. 

In Fig. 12 we show some current densities derived from the same data we used to create 
Figs. 11. In the left panel we show the total right chiral number density = (tP^PriP) = 
YlhihR- Using the cQPA correlation function we can write = Tr[5^Pa] in terms of 
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the on-shell distribution functions as follows: 



j°R,( k ) = J R 1 + hv k )n kh - (1 - hv k )n kh ] - hv k —$lf< h+ + — 9f/ 
2 L J Wk Wl- 



(6.9) 



where the vacuum contribution has been subtracted off. The result again displays the 
oscillatory behaviour driven by nonzero coherence distributions f c ± in near relativistic 
modes right after the transition. The right chirality of our chosen interaction term (6.2) 
leads to a temporal charge separation between fields ip and q after the transition. This 
effect produces the effective chemical potential for the field q mentioned above, and a net 
particle-antiparticle number density asymmetry j° = (■ijj^ ip) = Ylh^h = Ylh ^M^/f] f° r 
the field ip. With the vacuum part again subtracted off we find: 



We show the net asymmetry in the right panel of Fig. 12. Of course the total, combined 
fermion number of tp and q fields is conserved, because no net fermion number production 
takes place in the absence of sphaleron processes. 

6.1 Error analysis 

In the cQPA scheme the collision terms are computed to the zeroth order of the gradient 
expansion in time derivatives of the mass function m(t). The effects of the neglected 
gradient corrections are twofold: first, the gradient corrections to the constraint equations 
would broaden the singular phase space structure of the 2-point functions S K> [25]. As 
a result e.g. the self-energy functions in Eq. (2.49) (or in expressions (6.8)) would not be 
evaluated at exactly on shell, but rather convolved with distributions which are peaked but 
nonsingular functions at the on-shell momenta. However if, as is generally the case, the 
self-energy functions are smooth in the scale set by the width of the distribution functions, 
the correction from this effect should be very small. 

Second, and more importantly, there are explicit gradient corrections to the collision 
term involving the derivatives dfS and <9f£, which (after performing the zeroth order re- 
summation) are proportional to mass gradients. It is not easy to estimate the precise 
numerical error coming from neglecting these gradients. However, a naive estimate for the 
relative error in the collision terms is given by dimensionless gradient expansion parameter: 




Note that e depends both on time and on the momenta involved. In the example of this 
section the bulk asymmetry is generated at momentum |k| ~ 0.5 moo (see Fig. 12), and for 
these modes, in the middle of the kink profile (6.1) at t = 0, where the gradients are largest, 
we find e « 0.8. However, the error rapidly gets smaller as one moves away from the wall 
and it dies exponentially when |Ai| > r. Thus e alone does not provide a direct estimate 
of the error on the final results. Now, the quantities proportional to particle-antiparticle 
asymmetries are the most sensitive to variation of collision terms. In these cases the true 
error should be bounded by the product of e with the wall width and T: 6^ el ~ err. In the 



i°( k ) = n kh ~ n kh . 



(6.10) 
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Figure 13: Left-chiral quark chemical potential S yqL s=s \ ^r- as a function of time. All components 
of the quark field are assumed to be in kinetic equilibrium at all times. 

current example we then find 5^ el ~ 0.1, suggesting our results may be accurate to within 
10 per cent. For the "bulk" -quantities (not proportional to asymmetry) the collision term 
is not dominating the evolution and consequently we should expect a much smaller induced 
error in these quantities. To test these arguments we redid our calculations by adding an 
arbitrary correction to all self-energies of the form 

where a = m,c. Numerical results verify the analytical error estimates given above: the 
chiral number densities and charge densities shown in figure 12 and the chemical potential 
for the g-species shown in figure 13 were seen to change by the expected 10%, while in 
the excess of the total particle number density and the total particle antiparticle coherence 
solutions of Fig. 11 the changes were at level < 1%. We believe that these evaluations 
provide a generous upper limit for the true errors, because in reality not all gradient 
corrections work coherently in the same direction. 

Finally, we wish to emphasize that these gradient corrections described above affect 
only to the evaluation of the collision term; the flow term and thus the free theory evolution 
is exact in this approximation. Indeed, the set of distribution functions f a form a complete 
reparametrization of the components of zeroth moment integral of S < , so that even in the 
region of large gradients where their interpretation as singular phase space distribution 
functions would break, they would still present the free theory evolution exactly. 

6.2 A toy model for coherent baryogenesis 

We can construct a very simple toy model of baryogenesis using the above results. Let us 
assume that the field q in our example corresponds to a standard model quark, while Y> is 
some new field with the complex time- varying mass m. We do not provide any details of the 
phase transition characteristics, except of noting that the order parameter field giving rise 
to ip mass need not be related to the field in the interaction (6.2). While not needed for 
a qualitative picture, such details would of course be necessary in a more serious modeling 
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effort. First, if q is a quark field, it would be held in kinetic equilibrium by the strong 
gauge interactions. We can assume that only one of the quarks (say the top quark) or 
any number of them are coupled to if) through interaction (6.2); variations of this type 
can be modeled by changing the strength of the coupling y. In addition, the strong chiral 
anomaly equilibrates the total left- and right chiral quark chemical potentials [3] , whereby 
the effective chemical potential sourcing the baryon number production is reduced by a 
factor of two: 

C-^- (6-13) 

Given £^ the total baryon asymmetry produced by the Electroweak anomaly can be com- 
puted from [3]: 

dn B 



where r sp h = n^a^T^ is the Chern-Simons number diffusion rate across the energy 
barrier which separates N- vacua of the SU(2) gauge theory, where K sp h = 20 ± 2 [26]. 
The second term describes sphaleron-induced relaxation of the baryon asymmetry in the 
symmetric phase: 

A ^ = PCS = 3 5^ M5i+ 5^ Mli, ( 6 - 15 ) 

i i 

where fics is the Chern-Simons- and /it/ the leptonic chemical potential and constant A is 
a model dependent quantity of order unity. We shall use a value A = 15/2 corresponding 
to the Standard Model particle content with flavour equilibrated chemical potentials apart 
from the right chiral leptons (see e.g. ref. [3]), but the actual value is not relevant. It is 
now easy to integrate the equation (6.14) for the total baryon number produced: 

3 f'°° 

n B = -r sph y dtQe~ kBt , (6.16) 



where 



kB = ^^£. (6-17) 



Now using rjB = ne/n 7 w 7tib/s and 7710 = lO 10 //^, we find: 

rp POO 

mo« 1800( ) / dtQ(t)e- h ^, (6.18) 

where t = tra^ and we took g*$ = 110 for the number of entropy degrees of freedom. The 
back-scattering term kg = k^/nioo ~ 1.24 x lO^T/moo is too small to affect our results 
as suggested above. For the results displayed in figures 11-12 we chose the parameters 
to optimize the display of the qualitative features of coherence solutions in the plots. We 
redid the numerical analysis for the baryogenesis model taking to = 0, t = 3moo, A9 = ir/2 
for the kink, and T = 5moo, m q = 0, = bm^ and y = 0.8 for the interactions. The 
resulting chemical potential £ qL (t) is shown in Fig. 13, and the baryon asymmetry given 
by Eq. (6.18) is 

7/10-7.5. (6.19) 
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Let us point out that this model is not realistic in the sense that a single fermion complex 
mass term breaks the CPT-symmetry. The model should then be thought of being a part 
of a more complex multi-field mixing scenario where the single-field treatment can be used 
to approximate diagonal mass eigenstates in the non-degenerate limit. Another obvious 
concern is the speed of transition giving rise to the V'-mass. We do not try to speculate on 
this issue any further, except note that perhaps such a "quench" could be effected through 
a coupling to another scalar field along the lines discussed in ref. [27]. 

7. Conclusions and outlook 

In this work we have reformulated and extended our recently introduced quantum kinetic 
theory for coherent, interacting fermionic and scalar fields in spatially homogeneous and 
isotropic systems. Our formalism is based on the coherent quasiparticle approximation 
(cQPA), where nonlocal coherence information is encoded in new spectral solutions at 
off-shell momenta. We have used the Schwinger-Keldysh formalism of non-equilibrium 
quantum field theory with pertinent approximations to derive the cQPA propagators in the 
spatially homogeneous background. In particular, we have now introduced the familiar 4- 
dimensional Dirac notation for fermions, which allows us to use the standard Dirac algebra 
in the reduction of the matrix elements. The Dirac-structure of the coherence propagator 
living at the shell ko = in particular is expressible as products of the usual mass-shell 
projectors onto positive and negative energies, or alternatively as direct products of positive 
and negative energy eigenspinors. The coherence propagators are always rapidly oscillating 
functions of time. 

We have shown that in the cQPA the Kadanoff-Baym equations for the 2-point func- 
tions in the Wigner representation reduce to a closed set of extended quantum Boltzmann 
equations (qBE) for the cQPA mass- and coherence-shell distribution functions. We ob- 
serve that the collision integrals in these qBE's need to be resummed to all orders in 
gradients because of the rapid oscillations of the coherence-shell functions. We have per- 
formed this resummation and derived a set of generalized momentum space Feynman rules 
for the theory, including effective 4-dimensional propagators endowed with the coherence, 
as well as an interaction vertex rule for a Yukawa theory. As a result of the resummation 
the vertex rule displays a local energy non-conservation in association with the coherence 
propagators, however in a way that never violates the energy conservation at the level of 
complete self-energy diagrams. 

We have given several examples of the diagrammatic calculations in our formalism, 
including 1-loop fermionic and scalar self-energies and a 2-loop fermionic self-energy in 
a Yukawa theory. We have also applied the formalism to construct a simple toy model 
for baryogenesis, where a coherent fermionic field is interacting with a thermal background 
scalar field and a Standard Model quark. In this model the CP-violation resulting from the 
C-breaking complex mass term of the ^-field and the chirality breaking decay interaction 
leads to a nonzero chemical potential for the left chiral quarks, which then acts as a seed 
asymmetry biasing the Chern-Simons diffusion rate to produce baryons via the Electroweak 
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anomaly. The final baryon asymmetry created via this mechanism can be roughly of the 
correct magnitude to explain the observed baryon excess in the universe. 

In this work we have considered only the case of a single fermionic and scalar field. 
A generalization of the current analysis to multiple coherently mixing fermion and scalar 
fields, including full flavour coherent propagators, flavoured quantum Boltzmann equations 
and a full set of flavour dependent Feynman rules will be presented in ref. [16]. The 
natural applications of the formalism with multiple mixing fields include for example the 
neutrino flavour oscillations in the early universe and the resonant leptogenesis [8], where 
the quantum coherence between the neutrino flavours may play an important role in the 
dynamics. For a realistic application to electroweak baryogenesis the formalism of this 
work needs to be generalized to the case of a stationary planar symmetry, where the 
coherence lives at the shell k z = 0. Based on our earlier work [11], we expect this case to 
be highly analogous, although technically somewhat more complicated than the spatially 
homogeneous formalism presented here. Also that formalism need to be developed for 
many mixing fields. However, despite the fact that most realistic applications of the cQPA 
formalism appear to need a generalization to mixing fields, most novel aspects related to 
the new coherence shells and resummation procedures necessary to define the appropriate 
Feynman rules should be qualitatively very similar to the techniques developed here. 
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